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Abstract 

Phase transitions (liquid-solid, solid-solid) triggered by temperature changes are studied in free nanosized clusters 
of TeFs (SFe) with different negative charges assigned to the fluorine atoms. Molecular dynamics simulations at 
constant energy show that the charge increase from gF=0.1e to ?F=0.25e shifts the melting temperature towards 
higher values and some of the metastable solid states disappear. The increased repulsive interaction maintains the 
order in molecular systems at higher temperatures. 
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1. Introduction 

Small free clusters of atoms or molecules ex- 
hibit solid-like or liquid-like properties that differ 
from the properties of their bulk counterpart. It 
has been realized that the cluster structure (micro- 
crystalline or amorphous) [1,2] depends on the pro- 
duction method. Clusters made of AFq molecules 
(A = S, Se, Te, U, W) have already been examined 
in some detail by several groups [2-6] to establish 
the existence of various structures and transforma- 
tions between them. Some of the microcrystalline 
states coexist dynamically [5,6] in a given tempera- 
ture interval and those which are only locally stable 
phases disappear when the cluster size increases as 
was confirmed both experimentally [1] and theoret- 
ically [6] . Those states correspond to a partial order 
of the molecular axes of symmetry. The system be- 
comes completely orientationally ordered at very 



low temperatures. The transition rate between the 
ordered and disordered states can be retrieved from 
the potential energy surface (PES) of the system, 
[7]. It has been shown in [6] that clusters of the 
same numbers of TeF$ or SF% molecules exhibit 
different dynamics despite the same symmetry of 
the molecules. The reason is that the topography 
of PES in the case of SFe clusters is shallower than 
that of TeFe clusters. 

In the present study we explore the changes of 
the cluster PES due to the changes of the charge 
distribution in a single molecule. Our hypothesis is 
that the molecular polarization is changed by us- 
ing different production methods. In order to find 
out the influence of the charge changes on orien- 
tational order-disorder phase transitions we have 
simulated the temperature behavior of molecular 
clusters with the help of a constant energy molec- 
ular dynamics. 
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2. Interaction Potential and Computational 
Procedure 

The main feature of the intermolecular interac- 
tion is the dependence on the mutual orientations 
of the molecules. There are experimental indica- 
tions that AF 6 molecules can be considered as rigid 
octahedra to a reasonable extent. A small negative 
charge should be assigned to the fluorine atoms to 
account for the chemical bond [3] . 

The intermolecular potential is presented as a 
sum of atom-atom interaction (fluorine-fluorine, 
tellurium-tellurium, fluorine-tellurium) to account 
for the orientational anisotropy: 
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where a, (3 denote either a fluorine or a tellurium 
(sulfur) atom; r"? is the distance between an a- 
atom in the i — th molecule and a /3-atom in the 
i — th molecule; n is the number of the molecules in 
the cluster. The parameters a a p and e Q( 3 have been 
fitted to the experimental diffraction results, [2] . 

The Coulomb term accounts for the small nega- 
tive charge qF assigned to the fluorine atoms and 
the positive charge carried by the central tellurium 
atom, q te — 6qF, which ensures a neutral molecule 
at distances much larger than the molecular size. 
Here we compare the temperature-driven transi- 
tions for the case of qF = O.le and qp = 0.25e, 
where e is the electron charge. 

The charge qF has been computed using LCAO 
(linear combinations of atomic orbitals) with pla- 
nar basic functions (qp = O.le) and Gaussian basic 
functions (qp = 0.25e). 

The potential, Eq.l, has been used to solve the 
classical equation of motion written in the Hamil- 
tonian form with the help of a constant-energy MD 
method. The velocity- Verlet [8] algorithm with a 
time step of 5/s has been implemented. This is a 



step optimized in [6] to satisfy the requirements for 
long MD runs necessary to detect phase changes in 
a computer experiment. The heating (cooling) of 
the system is performed by velocity rescaling and 
consequent thermalization [3] . 



3. Results and Conclusions 

We have investigated the thermal behavior of 
clusters consisting of 89 TeF§ (SFq) molecules in 
the temperature interval (50 -j- 140 K). The clus- 
ters of TeF§ molecules melt above 125K (qp = 
O.le) and 130K (qp = 0.25e) as is seen from the 
caloric curves in Fig.l, Fig. 2 and the Lindemann 
criterion, Fig. 3. The Lindemann coefficient 8 < 
0.08 corresponds to a solid-like phase (5 < 0.1 for 
bulk systems) [11]. Melting is a discontinuous tran- 
sition even in the case of nanosized clusters. The 
heating and cooling of clusters (Fig.l) demonstrate 
a hysteresis, which signals a discontinuous transi- 
tion. The charge increase shifts the melting point 
towards higher temperature: in comparison to a 
cluster with less charged atoms, Fig. 2. [ the larger 
charge, the more robust is the cluster ] . 

The freezing point shifts towards lower temper- 
atures. 

The hysteresis area is larger for the case of 
qF = 0.25e than for qp = O.le. One could specu- 
late about a "larger memory" in systems having a 
" larger charge" . Another important distinction be- 
tween less and more charged fluorine atoms is the 
structure adopted by the clusters below the freez- 
ing point. The distributions N(cos9) of the mutual 
molecular orientations with qp = O.le, Fig. 4, show 
that they transform step by step from liquid to 
a partially ordered solid (phase A) to an ordered 
solid (phase B). The clusters with qp = 0.25e, 
Fig. 5, transform directly from liquid to an ordered 
state. The radial distribution g(r) [10], Fig. 6, is 
a diagnostic of the lattice structure adopted by 
clusters. We compute g(r) for the molecular center 
of masses which is insensitive to the molecular ori- 
entations. Fig. 6 was obtained as follows: starting 
from a low-temperature configuration obtained 
from simulations of a cluster with qF = O.le, we 
change to qF = 0.25e and heat the cluster until it 



2 



melts. Then the cluster is cooled to a solid state. 
g(r) arc plotted at the same T = IIOK on the 
cooling and heating branch. The structures are ob- 
viously different. Typical configurations obtained 
by quenching [9] of the MD trajectory are shown 
in Fig.7 [(a) - for q F = O.le; (b) - for q F = 0.25e]. 
In both cases we obtain two orientations popu- 
lated by the molecules but the arrangements of 
"layers" in the clusters arc different. 

We conclude that the charge change influences 
both the melting temperature and structural 
transformations. Meta-stable solid states with par- 
tial order disappear for the case of larger charge. 
The volume of the cluster decreases when the 
charge increases and the cluster is more robust on 
heating. 

Our results show that the charge increase shifts 
the transitions temperature towards higher val- 
ues and some of the metastable states disappear. 
This confirms our understanding that the repul- 
sive interaction maintains the order in molecular 
systems. 

Acknowledgments 

NATO Grant (CLG SA(PST.CLG.976363)5437 
is acknowledged. The work has been partially sup- 
ported by the University of Sofia Scientific Fund 
(2001). 

References 

[1] E. Lefterova, I. Dragieva, V. Tchanev, D. Mehandjicv, 
M. Mikhov, Journal of Magn. and Magn. Mat. 140-144 
(1995) 457-458 

[2] L.S. Bartcll and Shimin Xu, J. Phys. Chem. 95 (1991) 
8939. 

[3] A. Proykova and R.S. Berry, Z. Phys. D 40 (1997) 215. 

[4] J.B. Maillct, A. Boutin, S. Buttefey, F. Calvo and A. 
H. Fuchs, J. Chem. Phys. 109 (1998) 329. 

[5] S. Pisov and A. Proykova, in Proc. of the Meeting 
of Young Researchers in Physics, Heron Press 2001; I. 
Daykov and A. Proykova, ibit.. 

[6] A. Proykova, S. Pisov and R. S. Berry, J. Chem. Phys. 
(accepted: AIPID 515140JCP) 

[7] K.D. Ball, R.S. Berry, R.E. Kunz, F.Y. Li, A. Proykova, 
D. Wales, Science Vol.271 877-1024 (1996) 963-965 

[8] L. Verlet, Phys. Rev. 159 (1967) 98. 

[9] F. H. Stillinger and T. A. Weber, Kinam 3 (1981) 159. 



[10] M. P. Allen and D. J. Tildcslcy, Computer Simulation 
of Liquids (Clarendon Press, Oxford, 1994) p. 184 

[11] J. Jellinek, Metal-Ligand Interactions (Kluwcr 
Academic Publishers, Printed in the Netherlands, 1996) 
p.325-360. 



95 100 105 110 115 120 125 130 135 140 145 

Temperature [Kj 

Fig. 1. Caloric curve for TeFe(89) cluster, with qp = 0.25e 
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Fig. 2. Caloric curve for TeFe(89) cluster, with qp = O.le 
and qp = 0.25e 
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Fig. 3. Lindemann criterion for Xe_Fe(89) cluster, with 
qp = O.le and qp = 0.25e 
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6. Radial distribution g(r) for Tei<6(89) cluster, with 
= 0.25e 
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Fig. 4. Orientation distribution N(cos8) for TeFe(89) clus- 
ter, with qp = O.le 




Fig. 5. Orientation distribution N(cos9) for TeFg(89) clus- 
ter, with qp = 0.25e 
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Fig. 7. (a) A quenched configuration of TeF§ (89) clusters 
with qp = O.le: a sequence of two-one-two layers with 
specifically oriented molecules is seen; (b) A quenched con- 
figuration of the same cluster with qp = 0.25e. The cluster 
crystallizes in another structure: three rows in one orien- 
tation and two rows in another orientation. 
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